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Hadron Mass Extraction from Lattice QCD 
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The extraction of quantities from lattice QCD calculations at realistic quark masses is of considerable impor- 
tance. Whilst physical quark masses are some way off, the recent advances in the calculation of hadron masses 
within full QCD now invite improved extrapolation methods. We show that, provided the correct chiral behaviour 
of QCD is respected in the extrapolation to realistic quark masses, one can indeed obtain a fairly reliable deter- 
mination of masses, the sigma commutator and the J parameter. We summarise these findings by presenting the 
nonanalytic behaviour of nucleon and rho masses in the standard Edinburgh plot. 
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1. INTRODUCTION 

There are well-known difficulties associated 
with making dynamical fermion lattice QCD cal- 
culations at light quark masses. There is the need 
however, to relate quantities calculated on the lat- 
tice with physical observables, hence results are 
required at physical quark masses. These two 
mutually exclusive restrictions on the field have 
motivated the necessity for extrapolation from 
the region in which calculations are able to be 
performed — that is, the region of unphysically 
heavy quarks — to lighter masses, including the 
physical quark masses. In this paper we discuss 
the construction and application of an extrapo- 
lation method for masses that respects the 
correct chiral behaviour of QCD and also allows 
the extraction of other quantities . This ap- 
proach is not limited to the case of masses in dy- 
namical fermion lattice QCD calculations. Other 
successes of this approach may be found, for ex- 
ample, in the extrapolation of baryon charge radii 
p , magnetic moments structure functions 
jep and quenched QCD data g. 
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1.1. Goldstone Boson Loops 

It is accepted that Goldstone Boson loops play 
an important role in all hadronic properties — 
their role is in one sense the basis of Chiral Pertur- 
bation Theory (xPT). Lattice QCD calculations, 
as an ab initio approach to calculating quantities 
in QCD, implicitly includes these loop contribu- 
tions. It has become clear recently, with calcu- 
lations appearing at lighter quark masses |p|,|lo|, 
that the nai've linear extrapolation methods are 
not reproducing the data. In particular in |ic| ] it 
was stated 

"The existence of curvature [at 
small quark masses] is observed, ne- 
cessitating a cubic Ansatz for extrap- 
olation to the chiral limit." 

The following section reviews how the inclusion 
of chiral physics allows reliable extrapolations of 
lattice QCD calculations Section 3 reports 

new results for the Edinburgh plot. 

2. EXTRAPOLATION METHODS 

In QCD chiral symmetry is dynamically bro- 
ken, and the pion is almost a Goldstone boson. It 
plays a significant role in the self-energy contribu- 
tions to the N and A, because of the strong cou- 
pling to the baryons. Chiral symmetry requires 
that, in the region where perturbations around 
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light quarks makes sense, the mass of the nucleon 
has the form 

TOAr(m^) = m^^^ + ami + 13ml- + 

jmllnm^ + ... , (1) 

where mj^\ a, P and 7 are functions of the strong 
coupHng constant. In particular the values of the 
coefficients of the non-analytic (in quark mass) 
terms — recall that oc mq — are known ex- 
actly from xPT. However it is only recent results 
from the lattice that have indicated any need of 
higher order terms beyond that of a linear extrap- 
olation in quark mass (or mj). 

2.1. Chirally Motivated Form 

An attempt at having a chirally motivated form 
for extrapolating masses has been 

mNirriTr) — mo + ami + (3ml , (2) 

where mp, a and /? are fit parameters. Naively 
this is a good choice. It reflects the known non- 
analyticity from xPT and still reproduces the lat- 
tice results. The problem with this method is as- 
sociated with the choice of (3. The value of the 
coefficient of the cubic term is known explicitly in 
xPT. So a functional form, motivated by chiral 
symmetry, should preserve the known value of f3. 
Optimising f3 via a best fit to existing lattice data 
provides —0.55 GeV~^. However, the result from 
xPT is —5.6 GeV~^. That the coefficient is so 
small is not surprising. The functional form at- 
tempts to reproduce the lattice data over a large 
range of ml , where the data is predominantly lin- 
ear — as can be seen in Fig. |^. However xPT is 
an expansion about the massless quark limit and 
would not be expected to be applicable (or even 
convergent) at such large quark masses. 

2.2. Current Calculation 

It has been found ||l|,|| that by retaining the 
contributions to the self-energy of the hadron 
mass that vary the most rapidly with m^r near the 
chiral limit, a successful extrapolation method 
may be formed. This methodology includes the 
most important non-analytic structure in the 
hadron mass near the chiral limit with exactly the 
correct coefficients. The pion mass dependence of 



the masses of the N, A and p are: 

mN = flo +a2m^ +a-jvAr7r(AAr,m^) 

+0-JVA7r(Ajv, m^) (3) 

rriA = &o + b2ml + ctaAtt ( Aa , ) 

+fTAA'7r(AA,?Tl^) (4) 

rrip = Co + C2m^ + crp^^(Ap,m^) 

(Ap,m^) (5) 

where ctabc indicates the contribution from the 
A ^ BC ^ A self-energy process. The expres- 
sions for these self-energy contributions for the N 
and A may be found in The two significant 
processes for the p are the p — > lot: and p — > rnr 
self-energies and they are presented in ||| . 

An additional level of detail explicitly included 
in these extrapolation methods is the inclusion 
of the decay channels (in the case of the A the 
process A Ntt) . This process is not included in 
other methods, and yet is a vitally important and 
physically based consideration. However, because 
of the finite nature of the lattice, decays are not 
always possible. The finite periodic volume of the 
lattice restricts the available momenta to discrete 
values 

Stttt-i, , Lii Li, , , 

ku = — r-^, with - -f < ri. < -f 6 

where and a are the lattice size and spacing 
in the p, direction, respectively. 

Figure |l| indicates the expected behaviour of 
the masses of the N , A and p using Eqs. (^), (^) 
and (|^), with the physical masses being 0.940, 
1.173, 0.713 GeV respectively.^ We also present 
an error analysis of the fitting for the particular 
case of the p meson in Fig. H The shaded region 
is bounded below by an increase of Icr in the x^ 
per degree of freedom of the fit, and above by a 
physical constraint in our approach. It is clear 
that whilst the central value of the extrapolation 
gives an acceptable value for the physical mass, 
the uncertainties are large. A Gedanken exper- 
iment performed in |^ suggests that a ten-fold 
increase in the number of configurations at the 
lowest pion mass data point [ml '--^ 0.1 GeV^) 

^The excellent agreement with the experimental mass of 
the nucleon is coincidental. 




Figure 1. Two- flavour, dynamical fermion lattice 
QCD data for the A, N and vector meson (p) 
mass data from UKQCD j|] (open circles) and 
CP-PACS |l§ (filled circles). The solid lines are 
the continuum limit, infinite volume predictions 
of Eqs. (^, (H) and (^). The squares (barely dis- 
cernable from the data) are the predicted masses 
on a lattice of the same dimensions as the data 
at that pion mass. 



would reduce the uncertainty in the extrapolated 
value to the 5% level. 

3. OTHER QUANTITIES 

The advantage of calculating the mass of the 
hadrons in the manner described above is that the 
form allows the direct extraction of other proper- 
ties of the hadron that depend upon the quark 
mass dependence of the hadron mass. 

3.1. The Sigma Commutator 

The sigma commutator is a direct source of in- 
formation about chiral symmetry breaking within 
QCD As such it is a quantity of considerable 
importance to extract from lattice QCD calcula- 
tions. The form of the commutator is 

CTjv fh{N\uu + dd\N) (7) 



Figure 2. Analysis of the lattice data for the vec- 
tor meson (p) mass calculated by CP-PACS jlOj 
as a function of m J . The shaded area is bounded 
below by a la error bar. The upper bound is 
limited by a physical constraint discussed in |^ . 



where m is the average mass of the up and down 
quarks. 

(Tjv is not directly accessible via experiment, 
however world data suggests a value of 45 ± 8 
MeV 0. Early attempts at evaluating Eq. (||) 
found results in the range 15 to 25 MeV, and the 
attention soon changed to evaluating the matrix 
element, Eq. directly. In quenched calcula- 
tions the results were in the 40-60 MeV range, 
but a two flavour dynamical fermion calculation 
by the SESAM collaboration Q found a value 
of 18 ± 5 MeV. The difficulties associated with 
these approaches are two-fold. Firstly, the scale 
independent quantity of ctat must be constructed 
from the renormalisation depended quantities m 
and {N\uu + dd\N) . Additionally there still is the 
need to extrapolate the quantities to the physical 
pion mass. 

Our recent work showed that provided the ex- 
trapolation method is under control the evalua- 
tion of ctat at TTiTT = 140 MeV, is a straightforward 
calculation. The important advantage of this ap- 
proach is that one need only work with renormal- 
isation group invariant quantities. 

We discussed previously how a chirally moti- 
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vated form, Eq. (||), will not reproduce the lat- 
tice data if the coefficient of the term is the 
value required by xPT. However we also showed 
that allowing this coefficient to be a fit param- 
eter results in a value that is wrong by almost 
an order of magnitude. This becomes even more 
significant in the case of the sigma commutator. 
The required derivative promotes this coefficient 
to greater significance and the sign of the terms 
acts to reduce the value of <7Ar. However this is 
not an issue with the extrapolation forms dis- 
cussed above. The sign and magnitude of the 
cubic term is exactly that predicted by xPT, but 
the effects are countered by higher order terms 
— resulting in a prediction for the value of cttv 
that included the correct chiral physics. We find 
1^ that the value of the sigma commutator is ap- 
proximately 45 MeV. 

3.2. The J Parameter 

This dimensionless parameter was proposed as 
a quantitative measure, independent of the need 
for extrapolation — an ideal lattice observable 
lUll . It has the form 



J 



dirii, 



dvii 



rriK* 



m, 



(9) 



(10) 



K 



which, by substituting the experimental mass val- 
ues, yields the value Q 

J = 0.48(2) , 

In Fig. 1^ we present the value of the J parameter 
as obtained from Eqs. (||) and (||). The detailed 
slope of the curve is parameter dependent, how- 
ever the presence of the cusp is model indepen- 
dent. The cusp is a result of the two pion cut in 
the rho spectral function and has been ignored in 
previous attempts at evaluating the J parameter. 
We find a value for the J parameter of 45(7) in 
good agreement with the experimental value. We 
note, however, that if the point of evaluation cor- 
responded to ~ 0.15 GeV^ the J parameter 
would have been around 50% larger. 

3.3. Edinburgh Plot 

The baryon and meson masses on the lattice 
are all determined modulo the lattice spacing - a 



0.4 




Figure 3. The solid curve is a plot of the value 
of the J-parameter as a function of obtained 
from Eq. (||) and the best fit to the lattice re- 
sults. The vertical dotted line shows the point at 
which the J-parameter is evaluated (rfiplm-^ = 
1.8). The horizontal line displays the experimen- 
tal value (0.48) plotted between the physical val- 



scale that must be determined from some piece of 
data external to the lattice. One method of re- 
moving this scale is by plotting a ratio of masses 
— the Edinburgh plot. In Fig. ^ we present a pre- 
diction for the infinite volume, continuum limit 
extrapolation of the lattice data previously pre- 
sented. The two points known explicitly are indi- 
cated by open stars on the plot. The first known 
point is ratio of the physical masses of the tt, p 
and N. The second point is the heavy quark limit, 
when the masses of the hadrons become propor- 
tional to the constituent quarks. The effect of the 
opening of the decay channel of the rho is visible 
at rriTr/mp = 0.5. The effects induced, and the 
expected behaviour on the finite sized lattice will 
be presented in a future work p5|] . 

4. SUMMARY 

The importance of including the correct chi- 
ral behaviour in extrapolation methods is becom- 
ing more important as dynamical lattice QCD re- 
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Figure 4. Edinburgh plot for CP-PACS (filled 
symbols) and UKQCD (open symbols) results. 
The stars represent the known limiting cases, at 
the physical and heavy quark limits respectively. 
The solid line is the infinite volume, continuum 
limit behaviour predicted by our functional forms 
for the extrapolation of the N and p masses. 



suits appear at lighter quark masses. The suc- 
cesses of the approach outlined above include not 
only predictions for the physical masses of the 
hadrons investigated, but other quantities suc- 
cessfully reproduced. These other successes in- 
clude the sigma commutator and the J parame- 
ter — both of which have been a thorn in the side 
of dynamical fermion calculations. It is through 
the inclusion of the dominant chiral physics, the 
recognition that decay channels are important, 
and the understanding of some of the finite size 
lattice artifacts that we have been able to suc- 
cessfully extrapolate the Edinburgh plot to the 
known physical limit. 
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